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ABSTRACT 

We present Giant Metrewave Radio Telescope (GMRT) 240 MHz observations of the nearby 
luminous FR I radio source 3C 270, in the group-central elliptical NGC 4261. Combining 
these data with reprocessed Very Large Array (VLA) 1.55 and 4.8 GHz observations, we 
produce spectral index maps that reveal a constant spectral index along the jets and a gradual 
steepening from the ends of the jets through the lobes towards the nucleus. A Jaffe & Perola 
(JP) model fitted to the integrated spectrum of the source gives an asymptotic low-frequency 
index of ainj = 0-53lo'o2’ while JP models fitted to the observed spectral index trend along 
the lobes allow us to estimate radiative ages of ^ 29 Myr and ~ 37 Myr for the west and east 
lobes respectively. Our age estimates are a factor of two lower than the 75-Myr upper limit 
derived from X-ray data ( O’Sullivan et al.|[2011| l. We find unlikely the scenario of an early 
supersonic phase in which the lobe expanded into the ISM at approximately Mach 6 (3500 
km s“^), and suggest that either the source underwent multiple AON outbursts with possible 
large changes in jet power, or possibly that the source age that we find is due to a backflow 
that transports young electrons from the jet tips through the lobes toward the nucleus relatively 
quickly. We calculate that in the lobes the energy ratio of non-radiating to radiating particles 
is ^ 4 — 24 indicating significant gas entrainment. If the lobes are in pressure balance with 
their surroundings, the total energy required to heat the entrained material is 10^® erg, ~40% 
of the total enthalpy of the lobes. 


Key words: galaxies: groups: general - ageing, radio galaxies: individual (3C 270, NGC 
4261) - galaxies: active 


1 INTRODUCTION 

The study of radio galaxies has become increasingly important as 
they appear to be natural laboratories uniquely suited to the study 
of one of the most energetic phenomena in the Universe: relativistic 
jets arising from the associated active galactic nuclei (AGN). The 
reason this phenomenon is so important is that radio-loud activity 
of the central AGN is now considered to be the most likely source 
of the heating mechanism injecting energy into the hot gaseous ha¬ 
los of groups and clusters of galaxies, thereby balancing the effects 
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of radiative cooling i Fabian|2012[ [McNamara & Nulsen|2007[|Pe-| 


[terson & Fabian|2006|. 


The evolution of radio galaxies is commonly accepted to pass 
through three phases (as described in, e.g., [Kraft et al.|2007l l. In the 
beginning, the inflation of the lobes surrounding the jets is highly 
supersonic as the lobes appear to be significantly over-pressured in 
relation to the ambient medium. After this early, in most sources 
short-lived, supersonic phase, the expansion of the lobes continues 
and their internal pressure drops, approaching equilibrium with the 
ambient gas. In the final phase jets shut down and the lobes even¬ 
tually become unobservable due to energy losses. Examples of the 
complex interactions between radio lobes and a surrounding intra- 
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cluster or intra-group medium (IGM) are relatively common in the 
local universe. Radio galaxies in groups and clusters are observed 
to drive shocks into the IGM ( |Nulsen et al.|2005||Fabian et al.|2006 [ 
[Forman et al.|2007{[David et al.|20lT) which in some cases are de¬ 
tectable as X-ray surface brightness discontinuities (e.g. in 3C 310, 


[Kraft et al.|[2m2l and PKS B2152-699, [ Worrall et al.[[20!2l l, and 
to inflate cavities (e.g. [Birzan et al.|2008 O’Sullivan et al.|2011| 
[Fabian|2012y 


Present knowledge of the physical processes acting behind 
this phenomenology is still rather poor. In order to understand how 
AGN feedback works, the combination of high-quality X-ray data 
and high-sensitivity radio observations is essential. X-ray observa¬ 
tions provide information on the physical properties of the envi¬ 
ronment surrounding the radio source, most notably the pressure of 
the IGM that the radio lobes are expanding into, and allow dynam¬ 
ical estimates of the expansion timescale and AGN power output. 
Multi-frequency radio data provide information on the cycle of jet 
activity, including independent estimates of the source expansion 
timescale and indications of the physical mechanisms of energy in¬ 
jection ( [Giacintucci et al.|2008[[20T2) . 


With this aim we have chosen to examine NGC 4261, a nearby 
galaxy which has an active galactic nucleus that contains a super- 


massive black hole of (4.9±1.0)x 10 “Mq (Ferrarese, Ford & Jaffe 


[1996[ l and hosts the low-power FR I radio source 3C 270. It is the 
brightest elliptical in a poor group ( [Forbes et al.|2006l > projected 
onto the outskirts of the nearby Virgo cluster. Its proximity means 
that its features can be observed with good sensitivity and linear 
resolution, allowing detailed study of the energetics of lobe infla¬ 
tion. 


3C 270 is composed of symmetric narrow kiloparsec-scale 
twin jets ( [Birkinshaw & Davies[[1985) which previous modelling 
studies indicate lie close to the plane of the sky (63° — 76°; [Finer] 
[et al.[[200T| [Laing & Bridle[[2014^ , and which flare and broaden 
before extending tens of kiloparsecs into the intragroup medium 
(IGM) jWorrall et al.|201^ . The jets emanate from a compact ra¬ 
dio core coincident with the optical nucleus of the galaxy ^Jones[ 
[ & Wehrle|1997[ l, where optical spectroscopy has revealed a low- 
ionization nuclear emission-line region ( [Ho, Filippenko & Sargent[ 
[1997[l. In addition, a nearly edge-on nuclear disk (100-pc-scale) 
rich in dust, molecular an d atomic gas jjaffe et al.|1993 1 , lies or¬ 
thogonal to the radio jets ( [Tremblay et al.|2009 1 . Laing & Bridle[ 
([2013^ examined the 1.4-4.8 GHz spectral index of the base of the 


jets and found values of ~0.65 at ~ 10 from the nucleus, jumping 
to 0.72 at ~ 14 and then declining to ~0.6-0.65 at 50 . The spec¬ 
tral index is defined ns Sv oc where Sv is the flux density at 
the frequency u. 


In the X-ray band, ROSAT images of the NGC 4261 group 
show diffuse thermal emission and a non-thermal nuclear compo¬ 
nent at low energies jWorrall & Birkinshaw|1994| l. Chandra obser¬ 
vations reveal X-ray jets in the inner few kiloparsecs of the radio 
jets ( jGliozzi, Sambruna & Brandt|2003[rZezas et al.[2005[ l and sug¬ 
gest a synchrotron origin for the jet X-ray emission. 


More recent analysis of the Chandra and XMM-Newlon data 
in [0’Sullivan et al.] ( [2011^ shows much structure in the gas in which 
the galaxy is embedded. The AGN jets have inflated two lobes, 
producing cavities in the IGM and building up rims of compressed 
hot gas which almost entirely enclose the lobes. The cavity rims 
consist of material that has been swept up and compressed by the 
expansion of the lobes, weakly heated by compression with no ev¬ 
idence of shock heating ( [O’Sullivan et al.|2011[ (. NGC 4261 hosts 
a small cool core with radius ~ 10 kpc and a temperature of ~ 0.6 


keV, with the IGM temperature in the group being ~ 1.6 keV 
( [Humphrey et al.[|2009^ . The structure of the core has been dis¬ 
turbed by the expansion of the radio lobes and is not spherically 
symmetric ( [O’Sullivan et al.|[2011[ l, with wedge-shaped surface- 
brightness decrements that suggest the jets have driven out the IGM 
(or galactic gas) from conical regions around their first few kilopar¬ 
secs ( [Worrall et al.|2010] . 

Evidence of anisotropy in the globular cluster distribution has 
been seen ( [Giordano et al.|2005| [Bonfini et al^[2012[|D’Abrusc^ 

[et al.|20T^ while deep optical imaging has revealed a weak tidal 

tail to the north-west and a tidal fan extending southeast from 
the galaxy ( [Tal et al.|2009) . These disturbed features suggest that 
NGC 4261 went through merging or tidal interactions with another 
galaxy within the past 1-2 Gyr. 

In this paper we present a detailed low-radio-frequency spec¬ 
tral study of the galaxy NGC 4261. Previous studies have left the 
question of the radiative age of the source unclear and with our 
study will try to resolve this problem. Our work is based on archival 
data from the Giant Metrewave Radio Telescope (GMRT) at 610 
MHz and 240 MHz and the Very Large Array (VLA) at 1.55 GHz 
and 4.8 GHz. The primary goals of this work are to determine the 
radiative age of 3C 270 and to understand the energetics and dy¬ 
namics of the interaction between the relativistic plasma of the ra¬ 
dio lobes and the X-ray emitting gas following [O’Sullivan et al.] 
( [2011^ . The paper is organised as follows. In Sections an d§ we 
describe the GMRT and VLA archival data used and the data re¬ 
duction. The radio images are presented in Sectionj^ In Sectionl^ 
we derive the physical properties of the radio source. In Section^ 
we discuss our results and consider the development of the radio 
source in the context of its environment. The summary and conclu¬ 
sions are finally presented in Section]^ 

We adopt a distance to NGC 4261 of 31.3 Mpc, which gives 
an angular scale of 1" — 0.151 kpc. This is consistent with the 
distance adopted by [Worrall et al. i2010i and negligibly different 
from the distance adopted by O’Sullivan et al.[ ( [2011[ ( who used an 
angular scale of 1” = 0.153 kpc for their X-ray analysis. When 
correcting for redshift we adopt a value of z=0.007378 ( [Cappel-[ 
[lari et al.|2011) . The radio spectral index a is defined according to 
Sv oc 1 '““, where Su is the flux density at the frequency v. 


2 GMRT OBSERVATIONS AND DATA REDUCTION 

NGC 4261 was observed using the GMRT at 240 MHz and 610 
MHz with total integration time of 4.5 hours in each band. At 610 
MHz the data were recorded using the lower and upper side bands 
for an observing bandwidth of 16 MHz each. For the 240 MHz data 
only one correlator was used to record the data in a total observ¬ 
ing bandwidth of 16 MHz where only half of the data were useful. 
Data for both frequencies were collected with 256 channels at a 
spectral resolution of 125 kHz per channel. The data were analysed 
and reduced using the NRAO Astronomical Image Processing Sys¬ 
tem (AIPS) package. We summarize the details of the observations 
in Table where we report the observing date, frequency, total 
observing bandwidth, total time on source, half-power beamwidth 
(HPBW) of the full array and rms level (1 ct) at full resolution. 

The data were edited to remove bad antennas and other obvi¬ 
ous defects. Initially the phase and flux calibrators were checked 
and the data with extreme phase differences were removed. This 
procedure was also applied to our source. The data were calibrated 
and the flux density scale was defined using the amplitude cali¬ 
brators 3C 286 and 3C 147. The task SETJY was used to set the 
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Table 1. Summary of GMRT radio observations 


Project Code 

Observation date 

Frequency 

(MHz) 

Bandwidth 

(MHz) 

Integration time 
(min) 

hpewH pa 

(full array,^ X ", °) 

rms 

(mJy bearn”^) 

15KMH01 

2009 Feb 15 

610 

32 

270 

7.32x4.77, 76.62 

1 

15KMH01 

2009 Feb 15 

240 

16 

270 

15.30x11.00, 67.17 

1.6 



Figure 1. GMRT 240 MHz radio contours of 3C 270, overlaid on the R 
band optical image from the Digitized Sky Survey (DSS). The contour lev¬ 
els (black) start at 4cr and rise by factors of 2. The Itr level in the radio 
image is ~ 1.6 mjy beam“^. Labels indicate the individual components of 
the radio source. Restoring beam: 15.30 X 11.00, RA.= 67.17°. 

flux density scale to “best VLA values (2010)”. The flux densities 
set for 3C 286 are 25.84±0.7 Jy for 240 MHz, 14.54±0.2 Jy for 
1.5 GHz and 7.46±0.05 Jy for 4.8 GHz, while for 3C 147 the flux 
density at 240 MHz is 48.65±1.46 Jy. The calibrator flux densities 
agree with|Perley & Butl^ ( |2013| l to within 1% and the |Scaife~^ 
|H ealdH2012| scale to within 3% except for 3C 147 at 240MHz, 
where the two scales disagree (60.95±2.26 Jy for |Scaife & Heald| 
|2012| at 240 MHz). The point source 1330-1-251 was used as phase 
calibrator for both 610 and 240 MHz. The bandpass was calibrated 
using observations of 3C 286, with channels that presented ampli¬ 
tude ‘spikes’ removed. At 610 MHz about 25% of the data were 
edited out. In order to increase signal-to-noise per channel during 
RFI removal while limiting the effects of bandwidth smearing, the 
channels were averaged into 6 of ~2 MHz each at 610 MHz, and 4 
of ~ 1.5 MHz each at 240 MHz. 

Due to the large field of view of the GMRT, the field was split 
into multiple facets for imaging. 56 facets were created at 610 MHz 
(with cellsize 1.5”) and 85 facets were created at 240 MHz (with 
cellsize 3”) so that in each facet the imaging field still approxi¬ 
mates a plane surface. When recombined for our final image these 
give a field of view of ~ 1.2° x 1.2° at 610 MHz and ~ 3° x 3° at 
240 MHz. After editing again the averaged data, repeated cycles of 
deconvolution were performed along with careful phase-only self¬ 
calibration in order to image the data. Most of the remaining noise 
in our final images arises from confusion and calibration uncer¬ 
tainties. At the lowest frequencies the phase errors are expected to 
originate from rapidly varying ionospheric delays. 

The final images were corrected for the primary beam pattern 
of the GMRT using the task PBCOR in AIPS. The rms noise level 
(la) achieved in the final full resolution images is 1 mJy at 610 
MHz, and 1.6 mJy at 240 MHz (see Table[^. 


[Chandra, Ray & Bhatnagar]j2004| > estimate the typical uncer¬ 
tainty on GMRT flux measurements to be ~5% at 610 MHz and 
~8% at 240 MHz, including contributions from the uncertainty in 
the measurement of calibrator flux and uncertainties introduced by 
the telescope hardware. These are the levels of uncertainty which 
we quote in Table In addition, they quote an elevation depen¬ 
dent error on the scale of 2-4%, arising from changes in telescope 
response between the calibrator and source positions. These errors 
dominate over the systematic errors in the flux density scale. 

Examination of the 610-MHz maps revealed ghost images 
containing up to ~10% of the flux of the main image to the N 
and S of the true radio structure, presumably because of residual 
delay and amplitude errors, though no attempt to improve these 
corrections was successful. The 610-MHz image is generally sim¬ 
ilar to images with comparable angular resolution from the VLA 
(Section]^, but since the ghost images introduce false structure, 
we choose to use the 610 MHz data only for an estimate of the to¬ 
tal flux density, and not for spectral analysis of regions within the 
source 

We imaged the phase and flux calibrators to check the integrity 
of the calibrations. While these were as expected, the flux densities 
derived for 8 background sources in the 3C 270 field were found 
to be systematically low by about 20% at 240 and 610 MHz as 
compared with expectations based on literature measurements. At 
240 MHz all source structures appeared normal, and it was consid¬ 
ered appropriate to increase all flux densities, including those for 
3C 270, by 20%. The same procedure was applied for the total flux 
density of 3C 270 at 610 MHz. We note that this discrepancy can¬ 
not be resolved by using the higher 240 MHz calibrator flux for 
3C 147 suggested by the |Scaife & Heald| ( |2012| l scale, since this 
only affects one of the two calibrators, and only at one frequency. 


3 VLA OBSERVATIONS AND DATA REDUCTION 

From the VLA public archive we acquired and reanalysed the ra¬ 
dio observations listed in Table The data were calibrated and 
imaged in AIPS. In order to achieve the best possible noise level 
and uv coverage for the 1.55 GHz and 4.86 GHz VLA archival 
data, several datasets were combined for mapping. Individual ob¬ 
servations were analysed separately until the stage where the source 
and calibrator uv data were separated. The SPLIT files were com¬ 
bined using DBCON and the analysis proceeded with phase-only 
self-calibration and imaging. For the 1.55 GHz data, as can be 
seen in Table[^ four different observations were combined together 
(AP0077, AV0088, AH0343 A and AH0343 B) in order to produce 
the final image with a total integration time of 61 minutes and an 
rms noise level of 0.8 mJy (Array B&C in Table [^. At 4.86 GHz 
only the D-array observation was used and analysed as it was suffi¬ 
ciently long and has a resolution comparable to the combined 1.55 
GHz image. The properties of the resulting maps are given in Ta¬ 
ble [H 

The task LGEOM was used to rotate the data where needed. 
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Table 2. Summary of archival VLA obsei'vations 


Project 

Observation date 

An'ay 

Frequency 

(GHz) 

Bandwidth 

(MHz) 

Integration time 
(min) 

AP0077 

1984 Apr 23 

c 

1.55 

2x50 

25 

AV0088 

1984 Jun 05 

c 

1.49 

2x50 

16 

AH0343 A 

1989 Mar 15 

B 

1.57 

2x50 

10 

AH0343 B 

1989 Jul 28 

BnC 

1.57 

2x50 

10 

AL0693 B 

2007 Mai-18 

D 

4.86 

2x50 

81 



40.0 35.0 12 : 19 : 30.0 25.0 20.0 15.0 10.0 


Right ascension 

Figure 2. VLA 1.55 GHz radio contours of 3C 270 overlaid on the radio 
emission. The 1 cr rms noise is ~ 0.8 mJy beam“^. Contour levels start at 
4cr and rise by factors of 2. Flux Density: 18.3 Jy, Restoring beam: 15.78 x 
10.69, RA.=-57.44°. 



40.0 35.0 12:19:30.0 25.0 20.0 15.0 10.0 


Right ascension 

Figure 3. VLA 4.86 GHz radio contours of 3C 270 overlaid on the radio 
emission. The I cr rms noise is ~ 0.4 mJy beam“^. Contour levels start at 
4 (t and rise by factors of 2. Flux Density: 8.3 Jy, Restoring beam: 19.01” X 
11.17”, RA.= -55.35°. 

to match the VLA and GMRT observations for spectral index map¬ 
ping. Finally the images were corrected for the primary beam re¬ 
sponse using the task PBCOR. Both images have spatial resolution 
similar to that of the GMRT 240 MFlz image. We again adopted the 
“best VLA values (2010)” flux density scale. 


4 RADIO IMAGES 

Figure [T] presents the GMRT 240 MHz radio image of 3C 270 at 
a resolution of 15.3 x 11.0 overlaid on the optical image from 
the Digitized Sky Survey. The contour levels start at Aa and rise 
by factors of 2. The most prominent components of the source, that 
is, the two symmetrical bright jets and their surrounding lobes, are 


labelled. The VLA images at 1.55 GHz and 4.86 GHz are shown 
in Figures|^and[^respectively. The jets are more prominent on the 
higher-frequency maps. 

The radio jets extend in opposite directions for several tens 
of kpc from the host galaxy of NGC 4261, having a similar mor¬ 
phology at all frequencies and at different restoring beams. The two 
bright jets are continuous and inflate two fairly round radio lobes 
that surround the jets. However, as expected, the prominence of the 
compact core increases with frequency, and it is clearly detected at 
1.55 and 4.86 GHz, but not at 240 MHz. 

There are two morphological differences between the two jets: 
first, the eastern jet bends to the north whereas the western one 
is straight but noticeably broadens and brightens ~l-2' from the 
nucleus; second, the western lobe is slightly less extended than the 
eastern one. These differences may be caused by interactions with 
their surroundings, but we see no evidence of a denser environment 
on either side in the X-ray. |0’ Sullivan eTaTj l lSoTT) note that the 
change in the direction of the eastern jet could be misleading, since 
we are viewing a projected image of a three-dimensional structure. 

Table summarises the properties of the source and its 
components including flux density, the spectral index in the 240 
MHz—1.55 GHz frequency range and the apparent linear size. The 
largest projected linear size of the source at all frequencies is ~80 
kpc, with a total angular extent of ~ 530” along the east/west axis. 
The jets have similar flux density and size, with the eastern jet be¬ 
ing about 2 kpc longer, while the western jet flux density is ~10% 
greater than the eastern one. As is evident from Figures [T] and 
at lower frequencies the lobe emission extends far enough to 
connect the lobes north and south of the nucleus, but at 4.8 GHz 
no such connection is observed. This is the first clear indication of 
the strong spectral gradients present in the radio emission of NGC 
4261. 


5 RADIO SPECTRAL ANALYSIS AND PHYSICAL 
PARAMETERS 

The archival GMRT and VLA radio data were combined with data 
from the literature in order to calculate the integrated radio spec¬ 
trum of 3C 270. The physical parameters of 3C 270 were then cal¬ 
culated using the spectral index distribution of the source. The anal¬ 
ysis described below was performed using the Synage-l-l- package 
( |Murgia|2000| l. 

5.1 Spectral analysis 

The integrated radio spectrum between 22 MHz and 33 GHz for 
3C 270 is shown in Figure]^ Tablelists the literature flux den¬ 
sities used in the spectrum, and the source of each measurement. 
Adding our VLA measurements to the literature data, we fit a (|Jaffe| 
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Table 3. [Summary of reduced VLA images 


An'ay 

Frequency 

Integration time 

HPBW, PA 

rms 


(GHz) 

(min) 

("x",°) 

(mJy bearn”^) 

B&C 

1.55 

61 

19.0x11.2,-55.35 

0.8 

D 

4.86 

81 

15.8x10.7,-57.44 

0.4 


Table 4. Radio properties of NGC 4261 


Components 

S'240 MHz 
(Jy) 

Si. 55 GHz 
(Jy) 

54.8 GHz 

(Jy) 

^240MHz 

^1.55GHz 

^1.55GHz 

“ 4 . 8 GHz 

projected size 
(kpc) 

Total.... 

48.50±3.21 

18.30±0.92 

8.30±0.42 

0.41±0.08 

1.01±0.07 

80 X 40 

Core. 

0.48±0.04 

0.36±0.02 

0.37±0.02 

0.16±0.10 

0.03±0.08 

- 

West jet.. 

10.60±0.85 

3.80±0.19 

1.38±0.20 

0.47±0.09 

0.90±0.15 

37 X 7 

East jet.. 

8.00±0.64 

2.80±0.14 

1.13±0.06 

0.56±0.09 

0.80±0.07 

43 X 8 



Figure 4. Radio spectr-um of NGC 4261 between 22 MHz and 33 GHz. The 
black points are literature data, (see Table|^, and the red points are the VLA 
4.86 GHz and 1.55 GHz, and the GMRT 240 MHz analysed in this paper. 
The best-fit curve is a Jaffe-Perola (JP) model applied to the data. The best 
fitting value of the spectral index at injection, ai„j 

is given, as well as the reduced value from the fit. 


|& Perola|1973| l (JP) model to the integrated spectrum. Our aim in 
fitting this model is firstly to determine whether the low-frequency 
data can be well modelled as a power-law, whose slope will ap¬ 
proximate the injection spectral index, and secondly to allow us 
to compare our GMRT flux density measurements with the well- 
calibrated data from other frequencies. The model assumes that 
the radiative timescale of the electrons is much longer than their 
timescale for continuous isotropisation. The resulting fit describes 
the low frequency data well, with an asymptoti c low-frequency in - 


2005 I. 


dex of ainj = O.SSIq qj, typical of FR I jets i Young et al. 

The slope steepens above ~5 GHz, with a break frequency of ~85 
GHz, though this is likely driven to a high value owing to the high 
luminosity of the relatively young plasma in the jets and core. We 
note that a continuous-injection (Cl) model gives an identical fit. 


The maximum spatial scale resolved in our GMRT and VLA 
images is large compared to the source, with the possible exception 
of the 4.86 GHz VLA observation (~72', 28', 19' and 6' at 240, 
610, 1550 and 4860 MHz respectively). However, the VLA 1.55 
and 4.86 GHz data points lie close to the fit and are both in good 
agreement with previous measurements at these frequencies, indi¬ 
cating that the VLA data capture the great majority of the flux from 
the source at both frequencies. 

The GMRT data points are in good agreement with the model 
fit, suggesting that the 20% correction we applied to the GMRT flux 
calibration is correct. Without this correction, the GMRT 240 MHz 
flux density would fall ~2.8 (j below the model prediction. Includ¬ 
ing systematic uncertainties on the flux scale and the elevation- 
dependent error on GMRT flux calibration would bring the mea¬ 
surement to within la of the model. However, we note that all the 
bright sources in the GMRT fields for which multi-wavelength data 
are available show systematically low flux densities; while the dif¬ 
ference may not be strongly significant at a single frequency for 
3C 270, the agreement across multiple sources and both bands sug¬ 
gests that the problem is real and significant. 

To provide a more nuanced picture of the spectral index struc¬ 
ture within 3C 270, we created spectral index maps of the source 
in two frequency ranges, 240 MHz-L55 GHz and 1.55-4.86 GHz. 
Each set of images was produced in IMAGR using the same uv 
range, Gaussian taper and cellsize. We used a 20 circular restor¬ 
ing beam (chosen based on the lowest resolution 1.55 GHz im¬ 
age), corrected for the primary beam, and aligned the images. Re¬ 
gions where the flux density in either band was <5a significant 
were removed. The resulting spectral index maps are shown in Fig¬ 
ures and Typical uncertainties on the spectral index in the 
lobes (excluding systematic uncertainties on the flux density scale) 
are ±0.06 for the 1.55-4.8 GHz map and ±0.05 for 240 MHz- 
1.55 GHz. Uncertainties in the spectral indices of the jets and core 
can be found in Tablej?) 

In Figure]^ the unresolved core of the radio source can be 
seen. It has a spechal index 0 : 240 ° ~ 0.16 ± 0.10, probably as a re¬ 
sult of self-absorption (as indicated by the lack of an obvious core 
in the 240 MHz image. Fig. [TJ since the spectral index is rather 
frequency-dependent (high curvature). Moving outwards the spec¬ 
tral index becomes steeper and a constant value of a ~ 0.47±0.09 
can be seen in the jets, while the lobes that dominate the total in¬ 
dex in the lower-frequency band appear with constant but slightly 
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Figure 5. Map of the spectral index distribution between 240 MHz and 1.55 GHz based on images with 20 HPBW. The typical uncertainty on spectral index 
in the lobes is ±0.05. The values within 1 beam of the sharp edges of the lobes are likely to be unreliable. 
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Figure 6. Map of the spectral index distribution between 1.55 GHz and 4.8 GHz based on images with 20 HPBW (indicated by the circle on the lower left). 
The typical uncertainty on spectral index in the lobes is ±0.06. The values within 1 beam of the sharp edges of the lobes are likely to be unreliable. 


steeper values (a ~ 0.60 ± 0.09). In Figurej^the core is obvious 
with a flat spectral index Q 1550 ~ 0.03 ± 0.08. The steepening 
of the spectral index within the lobes is much clearer in this im¬ 
age as a spectral index of ~ 0.9 — 1 can be seen above and below 
the beginning of the jets. In the jets, the spectral index at the point 
where they are surrounded by the lobes is ~ 0.5 — 0.6 and remains 
constant throughout. The spectrum in the lobes gradually steepens 
towards the nucleus, reaching values of ~1±0.06. In both Figures 
5 and 6 , the features in the spectral-index maps at the edges of the 
source, and particularly at the ends of the lobes, are likely to be 


spurious and induced by residual reconstruction errors in regions 
of steep intensity gradient. 

We performed a fit of the observed spectral index trend 
throtmh the lobes using the techniques of |Murgi^ p003| l. As Fig¬ 
ures and Ushow, changes much more through the lobes 

than 0240 °, indicating that the break frequency is >1.55 GHz over 
most of the area. The Synage-l-l- package only considers spectral 
indices between two frequencies, and we therefore use it to model 
the trend in afifo, while confirming that the 240 MHz data at each 
position do not conflict with the higher frequency data. The spectral 
index trend was determined by averaging the flux densities of each 
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Figure 7. 1.55 GHz image in grey scale. Overlaid are the 1.55 GHz contour 
levels starting at 3.2 mJy bearn”^ and rising by factor of 2. The cyan cir¬ 
cles indicate the lobe regions used for the source age estimation using the 
spectral index profile method. 



Distance along the lobe (arcsec) 

Figure 8. Spectral-index distribution along the west lobe of the source cal¬ 
culated from the 1.55 and 4.86 GHz frequencies, derived using the circular 
regions shown in Fig. |7] The solid line represents the best JP model fit. 
The values of cx-inj and i^hreak along with the reduced from the fit are 
reported here. 

n 

frequency in 40 diameter circular regions starting in the region 
likely to contain the youngest plasma, at the end of the jets, and 
then moving back along the paths shown in Figure]^ to the oldest 
plasma toward the nucleus. The circular regions were chosen to be 
larger than the beam size so that the measurements of the spectral 
index are independent. Figures[^and|^show the resulting spectral 
index measurements in the west and east lobes respectively. 

The trends in the west and east lobe were fitted using a 
JP model. Assuming a simple constant source expansion velocity 
(d oc t), and knowing that the break frequency of the emission from 
a passively ageing relativistic plasma Vbreak oc 1/age^, we expect 
to find Vbreak OC d~^, where d is the distance from the end of the 
jet. The best fit for the observed spectral trend in west and east lobe 
of the source can be seen as a solid line in Figures[^and[^ respec¬ 
tively. We note that the model provides a reasonable description of 
the spectral index trend in both lobes, suggesting that our assump¬ 
tion of constant expansion velocity is not unreasonable. The model 
provided similar values of ainj and Vbreak for both lobe regions. In 
the west lobe the calculated injection spectral index from the model 
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Figure 9. As for Figure[^but for the east lobe. 

Table 5. Integrated radio fluxes for NGC 4261 drawn from the literature 


V 

(MHz) 

Flux Density 

(Jy) 

Reference 

22.3 

176 ± 12 

1 

26.3 

162 ± 21 

2 

38 

144 ± 7.2 

3 

86 

81.5 ± 1.1 

4 

178 

56.5 ± 2.8 

3 

750 

27.15 ± 0.23 

5 

1400 

18.62 ± 0.93 

3 

2695 

12.79 ± 0.64 

3 

5000 

8.26 ±0.15 

6 

10700 

4.69 ± 0.28 

7 

23000 

2.6 ±0.1 

8 

33000 

2.1 ± 0.2 

8 


REFERENCES.-! l)|Roger et al.|(T96^; (2) |Viner & Ericks^ (TTTS); (3) 


|Kellermann et al.U1969J; (4) Artyukh et al.H1969(; (5) Pauliny-Toth et al. 

fi^; (6) 

Pauliny-Toth & KellermannH1968f; (7)|Kellermann & Pauliny- 

|TothM1973 

; (8) WMARIBennett C. L. et al.H2003| 


was Oinj ~ 0.56 with a break frequency of Vbreak ~ 12.4 GHz 
while in the east lobe we used a fixed injection spectral index of 
ctinj ~ 0.5 and got a Vbreak ~ 7.5 GHz (see also Table|^. The 
break frequencies measured change little if the injection indices in 
the two lobe heads are assumed to be the same, either 0.5 or 0.56, 
and we note that both values are similar to the value of 0.53 found 
from the integrated spectrum. The similarity between the two lobes 
indicates that they have similar plasma properties, and we therefore 
consider it reasonable to fix ai„j in the east lobe. 

5.2 Physical parameters of 3C 270 

To calculate the key physical parameters of 3C 270, we have to 
make certain assumptions. Firstly, the relativistic particle and mag¬ 
netic field energy densities were assumed to be uniformly dis¬ 
tributed over the volume occupied by the radio source and in ap¬ 
proximate energy equipartition. A low energy cut-off ymin (where 
7 is the electron Lorentz factor) was introduced in the energy dis- 
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Table 6. Results of the spectral index profile analysis. Values mai'ked * were 
fixed during fitting. 


Region 

"inj 

^break 

"^rad 

'^adv 



(GHz) 

(10® yrs) 


West lobe 

U.DD_q Q3 

12.4l« 

29tl 

0.004c 

East lobe 

0.50* 

7 c+O.S 


0.004c 


tribution of the radiating electrons rather than the fixed 10 MHz — 
100 GHz frequency interval used in the classic equipartition equa¬ 
tions ( |Pacholczyk|1970[ >. We assumed that the magnetic field is un¬ 
ordered along the line of sight and used 7min = 100 in the electron 
energy distrihution, which corresponds to ~ 50 MeV. 

We used the injection spectral index derived from the west 
lohe ai„j=0.5, since this is in hroad agreement with the value for 
the east lohe and the source as a whole. The flux density at 240 
MHz (Table was used to calculate the radio luminosity of the 
source, since among the frequencies used in this work it is least 
affected by electron spectral ageing. As the lobes are not well de¬ 
scribed by simple ellipsoids, we divided them into a number of 
rectangular regions chosen to approximate the breadth of the lohe 
at a range of distances along the jet axis. These are shown in Fig- 
ure|10| Excluding the jets from the equipartition calculations for the 
lobes, we assume rotational symmetry about the jet axis to estimate 
the volume in each cylindrical region. We then group the regions, 
numbering these groups 1-8, to match the radial binning of the IGM 
pressure profile derived from X-ray observations. We also assume 
that the axis of the radio galaxy is in the plane of the sky. If the jet 
axis is oriented at 76° to the plane of the sky (as suggested hy the 
most recent modelling, |Laing & Bridle|2014| >, our estimated vol¬ 
ume will increase hy ~10%, altering our estimate of equipartition 
magnetic field by ~10%. 

The volumes, flux densities, spectral indices and equipartition 
magnetic field in regions 1-8 are listed in Table We also esti¬ 
mated the pressure of the relativistic plasma in each region, which 
is defined to be 


He, (1 -b k)E, 


2/roe 




( 1 ) 


where is the energy of the electron population, V is the 
volume, /To is the permeability of free space (47r in cgs units or 
47r X 10“^ in SI), (f) is the filling factor of the plasma (assumed 
to be 1) and k is the energy ratio of the non-radiating particles 
to the radiating relativistic electrons. We assume fc=l, indicating 
equal energy in radiating and non-radiating particles. The factor e 
describes the ordering of the magnetic field and we adopt e = 3, 
indicating a tangled magnetic field, as in |0’Sullivan et al^j2011^ . 

We find Bf.q ~5 /rG for both eastern and western lobes, 
though the outer parts of the lobes (regions 4 and 8) have higher 
Be,, close to the estimate for the source as a whole ~6.4 /rG. For 
this field (see later), electrons at 7^i„=100 radiate at ~5 MHz. 


5.3 Radiative age 

When considering the development of the radio jets and lobes, 
we adopt a simple conceptual model, in which at any given time 
plasma is rapidly transported out from the central engine to the tip 
of the jet. At this point the jet flow is halted by contact with the 
IGM, and the plasma driven out of the jet moves out radially to 
form the lobe. In this model, the age of the lobe plasma is linked 



Figure 10. Regions used for the estimation of the volume and the other 
physical parameters described in Table |7] The 3.2 /rJy (4a) contour is 
shown to indicate the limit of the detected lobe emission. 


to the distance along the jet; the back of the lobes (close to the nu¬ 
cleus) formed when the jet only extended a few kiloparsecs, while 
the tips of the lobes have only formed recently. In principle, back- 
flows could develop, transporting younger plasma away from the 
jet tips into regions of the lobe formed at earlier times. However, 
in the absence of such a flow, or external forces causing reaccel¬ 
eration, plasma in the lobes is expected to age passively, with the 
high-frequency break in its spectral index falling to lower frequen¬ 
cies over time. 

The radiative age of the radio source of NGC 4261 can be cal¬ 
culated if the break frequency of the source is known (e.g., |Myers| 
|& Spangler|1984| ). The magnetic field is considered to be uniform 
and constant over the whole of the source lifetime when radiative 
losses dominate over expansion losses. Neglecting reacceleration 
processes, the radiative age (the time since the last period of injec¬ 
tion of relativistic electrons began) can be obtained as 

^ 0.5 

frad = 1590 -[(1 + Z)l^breakl“° ® Myr, (2) 

where inbreak is expressed in GHz, and Be, and Bomb in ^G 
( [Parma et al.||2007) . Beq is the equipartition magnetic field and 
Bcmb = 3.2 is the equivalent magnetic field strength of the cos¬ 
mic microwave background (CMB) radiation at zero redshift, and 
equation thus includes both inverse Compton and synchrotron 
losses. 

Using equation]^ and the mean magnetic field for the source, 
5 /iG, we estimated a radiative age of ~29 Myr for the western 
lobe and ~37 Myr for the eastern one. We also derived an estimate 
of the advance speed of the lobe tip, v^d^/c, where c is the speed 
of light. This velocity was calculated as t)adv= LD/rrad, assuming 
a simple constant velocity and using the linear distance (LD) from 
the back of the lobe to the tip of the jet as measured from the radio 
images in Section]^ We found that for both lobes Vadv ~ 0.004c 
(~1300 km s“^ for the western lobe and ~1200 km s“^ for the 
eastern one). 

By comparison, the sound speed in the IGM in the environ¬ 
ment around NGC 4261, estimated from the measured X-ray tem¬ 
perature, increases from about 380 km s“^ near the core to 590 
km s“^ beyond ~100". The advance speed of the lobe tip there¬ 
fore appears to be supersonic relative to the IGM. Conversely, the 
Alfven speed in the lobe is about 10® km s“^ (using the lobe plasma 
properties deduced in Section [6^ . That is, speed of expansion of 
the lobe is strongly sub-Alfvenic, and so will not excite magnetic 
shocks in the lobes. The internal characteristic speed of the waves 
in the magnetized relativistic plasma were calculated from 
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Table 7. Physical parameters of 3C 270 and its components 


Region 

(arcsec) 

Tmax 

(arcsec) 

Vtot 

(kpc3) 

S 240 MHZ 

(mJy) 

>Sl550MHz 

(mJy) 

*^obs 

L 1550 

(10^2 W/Hz) 

Beq 

(pG) 

-^radio 

(10-^^dyn cm-2) 

West Lobe 

1 

10 

33 

1724 

998 ± 80 

316 ± 16 

0.62 

3.6 

5.0 

1.42 

2 

33 

63 

3861 

2067 ± 165 

650 ± 33 

0.62 

7.3 

4.9 

1.36 

3 

63 

168 

12081 

8741 ± 700 

3138 ± 157 

0.55 

35 

5.5 

1.73 

4 

168 

242 

2551 

2262 ± 181 

897 ± 45 

0.50 

10 

6.0 

2.06 

East Lobe 

5 

10 

33 

1056 

408 ± 33 

146 ± 7 

0.55 

1.6 

4.6 

1.21 

6 

33 

63 

1646 

976 ± 78 

287 ± 14 

0.66 

3.2 

4.9 

1.39 

7 

63 

168 

9753 

5868 ± 469 

2310 ± 116 

0.50 

26 

5.3 

1.64 

8 

168 

282 

4668 

5755 ± 460 

2202 ± 110 

0.52 

25 

6.5 

2.43 


- - : .. - ... 

I _j_ radio 

Y "^Pmag 

where va is the Alfven speed expressed in m s“^, Uptci is the en¬ 
ergy density of the relativistic particle population, and Pmag is the 
magnetic field pressure. 


6 DISCUSSION 

6.1 The age of the source 

Before considering the implications of our radiative age measure¬ 
ments, we first discuss potential sources of error in our approach. 
We have performed a relatively simple modelling of the 1.55- 
4.86 GHz spectral index trend in the lobes, and our model rests 
on a number of assumptions. These include constant source expan¬ 
sion velocity, constant magnetic field, our choice of spectral model, 
and that reacceleration, adiabatic expansion and mixing are not im¬ 
portant processes in the lobes. As mentioned in Section o the 
fact that our model provides a reasonable description of the spec¬ 
tral index trend suggests that the assumption of a constant expan¬ 
sion velocity is at least reasonable; we see no evidence that velocity 
changes are required, though of course other factors could disguise 
deviations from the trend. As a further test, we refitted the spec¬ 
tral index profiles, adding one additional region closer to the AGN. 
In both lobes, the fits were poorer, primarily because of the differ¬ 
ences in the extent of the emission seen in the lobes close to the 
nucleus (back of the lobes) at 1.55 and 4.8 GHz, which lead to di¬ 
vergence from the smooth trend in spectral index. For example, in 
the west lobe, the added region is close to the edge of the lobe at 
4.8 GHz, leading to a steeper than expected spectral index and a 
reduced x^=8.7. Both fits find lower break frequencies ~4-5 GHz, 
suggesting radiative ages ~45 Myr, older than our best estimate, 
but not by a large factor. 

Our use of the JP model may be over-simple, but we are lim¬ 
ited in this by the available data. |Harwood et aT1j2013^ compare the 
utility of several spectral models to describe spectral ageing in FR- 
II galaxies and conclude that the |TribbI^jl993| > model provides the 
best combination of fit quality and physical realism; Harwood et 
al. consider the JP model at least as physically accurate, but it pro¬ 
vides a poorer fit. However, they find that to differentiate between 
models requires data at 7-9 frequencies, whereas we have only 3 
available. Further observations are therefore needed to determine 
whether a different model might provide a better description (and 
age estimate) for 3C 270. 


Losses caused by adiabatic expansion would affect our age 
estimate, lowering the break frequency and therefore the esti¬ 
mated age. However, adiabatic losses decrease the luminosity of 
the plasma, and are more effective in doing so than they are in re¬ 
ducing Vbreak- We might expect losses in the jets as they flare from 
the narrowly collimated structures seen in the X-ray to the broad 
jets we observe in the radio, but the fact that the lobes have high 
surface brightnesses and smooth surface brightness distributions 
suggests that adiabatic losses within the lobes are probably not im¬ 
portant. Variation in the magnetic field within the lobes could affect 
our results, but on large scales we only measure a relatively small 
(~20%) change in field strength across the lobes (see Table]^, per¬ 
haps as little as 8% in the radius range used to estimate age. This 
is insufficient to balance the observed spectral index trend, and we 
therefore consider our assumption of a constant field to be a rea¬ 
sonable approximation. Small-scale field variations cannot be ruled 
out, but would again require improved spectral modelling and ad¬ 
ditional observations. 

Moving on to the results, our spectral-index profile modelling 
(Figures and 1^, implies ages of trad ~ 29 and 37 Myr for the 
west and east lobes respectively (Tabl e[^. These ages a re less than 
the <75 Myr upper limit estimated by |0’Sullivan et ah] ( |201 1| ) from 
the IGM sound speed. While the upper limit and our values are 
consistent, the factor two difference between the values seems to 
be significant. |Eilek| ( [l996^ first noted that the synchrotron ages for 
FRI and FRII radio sources derived from the spectra differed from 
the dynamical ages as much as a factor of 10. By this standard 
the factor of 2 discrepancy we observe is relatively mild. Several 
studies suggest plausible explanations for the age difference be¬ 
tween spectral and dynamical ages ([Blundell & Rawlings] 2000l>, 


including in situ acceleration of particles (e.g., Eilek|1996 Carillij 


jet al.|199T l and magnetic fields that vary with position in the lobes 
(e.g., Katz-Stone et al.|1993[|Eilek & Arendt|1996| l. Studies have 
shown much better agreement between dynamical and spectral ages 
in young sources (<10 Myr,[Matthews & Scheuer|1990||Kaiser etj 
[al.|1997[ , but the cause of the discrepancy in older sources is still 
unknown. We will discuss three possible explanations for the age 
discrepancy in 3C 270: i) a period of highly supersonic expansion 
early in the AGN outburst; ii) multiple outbursts or large changes 
in jet power; iii) a backflow capable of mixing younger electrons 
from the jets into the old plasma at the back of the lobes. 

In scenario (i), the expansion of the lobes was much faster 
in the past, so that the assumption of a constant near-sonic expan¬ 
sion velocity by [0’Sullivan et aL]j2011^ leads to an overestimate of 
the age of the source. For the larger eastern lobe, the radiative age 
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(37 Myr) implies that the lobes should be driving into the IGM 
at Mach~2, if they have expanded at constant velocity. Shocks 
in the IGM, driven by this expansion, would be detectable in the 
available X-ray observations as increases in temperature and sur¬ 
face brightness, and are ruled out by the measurements described 
in |0’Sullivan et al!] ( |2011[ l. However, there is also the possibility of 
that a period of supersonic expansion occurred early in the life of 
the source, followed by a period of slower expansion which we are 
now observing. In this case, the shocks would have detached from 
the lobe tips as the expansion speed dropped, and moved outwards, 
potentially passing out of the field of view of the XMM-Newton ob¬ 
servation. Assuming that the shock front is just outside the XMM 
field of view (~14' or ~130 kpc) we can calculate the duration of 
the supersonic phase. The minimum Mach number of such a shock 
would be Admin ~6 (~3500 km s“^). The corresponding maxi¬ 
mum period of supersonic expansion would be ~7.5 Myr, during 
which time the lobes would expand to ~28.5 kpc along the jet axis. 
After this, expansion would continue at the sound speed for the re¬ 
maining 3/4 of the lifetime of the radio source. 

A Mach~6 shock in the IGM would cause an increase in X- 
ray surface brightness of a factor ~13.5 and a temperature increase 
of a factor ~12. Even if the efficiency of energy deposition is low, 
such a shock would have strongly heated the IGM, and the observed 
temperature would be the result of this heating. For a 10% heating 
efficiency, we would expect the shock to have deposited ~10®® erg 
in the central 10 kpc alone, and the current rate of radiative energy 
loss, ~10"^^ erg s“^ would only have removed ~1% of this energy. 
This implies a pre-shock core temperature as low as 0.1 keV, lower 
than is observed in any group or cluster, and more comparable with 
the temperature of the interstellar gas in our own Galaxy than that 
expected in a giant elliptical. Since NGC 4261 appears to be fairly 
typical of the population of group and cluster-central FR-I radio 
galaxies, we might also expect to see a population of such sources 
driving shocks of similar strength, with evidence of strong shock 
heating. No such population has been identified. The scenario of a 
single, highly supersonic outburst therefore seems unlikely. 

In scenario (ii) the AGN would have undergone more than one 
outburst, or the jet power may have varied significantly over time. 
In this case our estimated age could correspond only to the newest 
outburst or period of high jet power. A restarted jet, or period of 
more powerful outflow, could drive shocks into the lobes, causing 
in situ reacceleration in the old plasma while producing new ex¬ 
pansion at the lobe tips. |Worrall et al.| ( |2010^ note structure that 
could be caused by variations on ^10"* yr timescales in the X-ray 
jet, while the changes in radio brightness along the western jet sug¬ 
gest variation on longer timescales. This scenario is plausible, and 
difficult to rule out with current data. 

In scenario (iii) a backflow would transport young, energetic 
electrons back towards the nucleus, reducing the mean age of the 
plasma. Numerical simulations (e.g., |Perucho & Marti|2007) ) sup¬ 
port the possibility of mildly relativistic backflows in a sheath sur¬ 
rounding light FRI jets. Modelling of the base of the jets in 3C 270 
does not reveal any evidence of such a backflow jLaing & Bri-| 
|dle|2 014), but it cannot be ruled out on larger scales. However, if 
backflows are limited to a sheath around the jets jLaing & Bridle| 
|2012| >, the question of how electrons are mixed outward into the 
larger lobe remains. A bulk backflow, with transport of and mixing 
of plasma throughout the lobe would be required, but seems less 
plausible given the lobe morphology. 

Considering the merits of these three possibilities, we favour 
scenario (ii), and suggest a model for the expansion history of the 
source in which, 75 Myr ago, an outburst swept its way through 


the X-ray gas, beginning the excavation of the cavities. These cav¬ 
ities grew in stages as the jet power varied over time, with the lobe 
plasma being refreshed by each new outburst or period of high jet 
power, the most recent of which is only ~30 Myr old. The overall 
lobe/cavity inflation time would thus be similar to the X-ray esti¬ 
mate, but the plasma we see in the lobes is a mixture of old and 
new particles injected by a period of stronger jet activity over the 
last ~30-40 Myr. 

6.2 Pressure balance and particle content 

The non-thermal energy density in the lobes of an extended ra¬ 
dio source cannot be determined by the radio data alone. Only a 
minimum value can be estimated when we assume the contribution 
from the magnetic field and the relativistic particles equal ( |Bur- | 
|bidge|[T958^ . The additional information needed to constrain the 
internal energy density can be obtained either from measurement 
of the inverse-Compton emission caused by the up-scattering of 
photons off the relativistic electrons in the lobes, or by assuming 
pressure equilibrium between the lobes and their environment and 
measuring the IGM pressure using X-ray data jLongairetal.|197^ . 

We first test our ability to detect inverse-Compton emission 
from the lobes. Assuming equipartition, |0’Sullivan et al.| j2011| | 
estimated an expected X-ray flux density from up-scattered cosmic 
microwave background (CMB) photons to be ~2 niy from each 
lobe at 1 keV. Using our improved constraints on the properties of 
the lobe plasma, we repeat this calculation, and find that we expect 
a flux of ~6.5xl0“'^'‘ erg s“^ in the 0.5-7 keV band (~4.1 njy 
at 1 keV) from each lobe. We then reprocessed the longest XMM- 
Newton dataset available (ObsID 0502120101, ~70 ks after clean¬ 
ing) using the latest calibration and most recent version of the XMM 
Science Analysis System (SAS 13.5). Following the techniques 
described in O’Sullivan et al., we find no detection of inverse- 
Compton flux from the lobes, but are able to place 3 (j upper lim¬ 
its on the 0.5-7 keV flux of <6x10“^'^ erg cm~^ s“^ (<3.8 nJy) 
for the west lobe, and <6.2x 10“^"^ erg cm~^ s“^ (<4.0 nJy) for 
the east lobe. These limits are marginally below the predicted flux, 
and so may suggest that the lobe plasma is out of equipartition, 
with more energy in the magnetic field than in the particle popula¬ 
tion. However, the uncertainties in the radio fluxes, XMM calibra¬ 
tion (see, e.g. |Nevalainen et al.|2010] |Tsujimoto et al.|2011^ and 
the choice of X-ray spectral model only weakly constrain depar¬ 
tures from equipartition in the sense of domination by the electron 
(and positron) energy density. Deeper X-ray data are needed to pro¬ 
vide a lower flux limit and determine whether the lobes are truly in 
equipartition. 

We also estimated the likely optical flux in the inner part of 
the lobes from up-scattering of infra-red photons emitted by the 
stellar population. Unfortunately, the expected signal is low, a fac¬ 
tor ~10“^ below the optical luminosity of the stellar population 
in the same region. We are therefore unable to constrain the prop¬ 
erties of the lobe emission from its inverse-Compton emission, at 
least with the currently available data. 

As inverse-Compton emission does not provide a definitive 
constraint, we compare the apparent lobe pressure with that of the 
surrounding IGM. Figure[TT]shows the pressure profiles of the ther¬ 
mal gas (Pgaa) along with the profile of the radio pressure derived 
in this study, as a function of the distance from the source’s cen¬ 
tre. The X-ray pressure profile is drawn from the XMM-Newton 
analysis of |0’ Sullivan et ah] j2011^ . As can be seen from Fig. [m 
the X-ray pressure profile declines with radius as expected, while 
the pressure profile derived from the radio lobes rises with radius. 
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Figure 11. Deprojected pressure profile of X-ray thermal pressure and radio 
pressure with radius from the center of the source. 


The radio lobes appear to be underpressured with respect to the sur¬ 
rounding environment in all regions and in both lobes by a factor 
ranging from ~2 to ~30. We note that the |0’Sullivan et 1^ 

IGM pressure profile is in excellent agreement with the earlier mea¬ 
surements of |Croston et al.| ( [2008| l, but that Croston et al. found a 
mean internal pressure for the radio lobes of ~9x 10“'^® erg cm“®, 
slightly lower than the pressures we measure in the lobes. 

A similar pressure imbalance has been found in previous lobe 
studies of other FR I radio sources, where Pgas/Pradio is usually 
>1, reaching values up to 100, assuming the standard equipartition 
arguments (e.g. Feretti et al.|1992[ I3unn et al.|2005[[Croston et al.| 


|2008[ [Birzan et al.|2008 Dunn et al.|2010| i. This can only suggest 

two possible situations for the radio lobes in the case of 3C 270: 1) 
Deviation from equipartition within the radio lobes with a higher 
magnetic field energy density than calculated, or 2) the equipar¬ 
tition argument is still valid but either the radio lobes contain a 
significant amount of non-radiating particles that contribute to the 
internal pressure or the filling factor cj) of the plasma is less than 
unity. Since we believe that jets entrain material in FR I galaxies 
and in this specific galaxy we have reason to think the filling factor 
might be ^ 1 at the back of the lobes (volume overestimation in¬ 
cluding part of the brighter X-ray emitting gas near the cool core), 
we examine the implications of the lack of pressure balance in the 
case where equipartition stands. 

Assuming that the radio lobes of 3C 270 are in approximate 
pressure equilibrium with the IGM, we have the ability to constrain 
the particle content of the lobes from the comparison of the appar¬ 
ent pressures of the relativistic plasma of the radio lobes and that of 
the thermal plasma of the IGM (e.g.,| Hardcastle & Worrall|2000 [ 
[Birzan et al.|2ro8| l. 

In Table the results shown for [1 + k)/(j) estimation were 
calculated by forcing Pradio to equal Pxray (note that in Table 
we adopted a fixed (1 + k)/(j> = 2). The regions of the radio pres¬ 
sure in this study were chosen so as to match the X-ray annuli used 
by|0’Sullivan et ^(|2011). Both lobes show the same trend with 


the (1 + k)/(f> factor rising as we move from the end of the jets 
(regions 4 and 8) to the back of the lobes near the core (regions 1 
and 5). The size of the bins were not exactly equal so the increase 
of(l + fc)/(/!)is not linear. 

Assuming that the filling factor 0=1, the additional pres¬ 
sure required by the rise of (1 + fc)/0 in regions 2-4 and 6-8 can 
be provided by non-radiative particles in the relativistic plasma of 
the lobes. This non-radiating component could be introduced via 
the entrainment and heating of thermal plasma by the radio jets. 
A probable source of entrained material is stellar mass loss from 
stars within the jets jPhinney||1983[ [Komissarov||1994[ |Hardcas-| 


|tle et al.|2003|[2007t . In |0’Sullivan et al.H2011t the rate of stel¬ 

lar mass-loss from AGB stars into each jet was estimated to be 
1.67x10“^ Mq yr~^. By assuming a 37 Myr timescale for the 
AGN outburst, we expect a minimum of ~60000 Mq of stellar 
material to be entrained by each jet. 

Thermal plasma entrained from stars within the radio jets rep¬ 
resents a plausible minimum mass of entrained material. Taking 
the limiting assumption that in regions 4 and 8 the only entrained 
material is that from stars in the jets (i.e., that no additional en¬ 
trainment occurs across the lobe/IGM boundary), we can estimate 
an upper limit on the temperature to which this thermal compo¬ 
nent must be heated in order to bring the lobes into thermal balance 
with their surroundings. We find the temperature of the entrained 
material to be ksT<l 1 MeV which corresponds to a plasma with 
electrons having a Lorentz factor 7<20, well into the relativistic 
regime. Electrons of 7~20 are expected to radiate in the equiparti¬ 
tion magnetic field at ~ 1 kHz, making them invisible in the observ¬ 
able radio band. Comparable values have been found by applying 
similar arguments to Cen A (|Wykes et al.| 2013l ) and 3C 31 ( |Croston| 
|& Hardcastle|2014l l. The very low density of this entrained thermal 
plasma component (n ~ 1.4 x 10“^ cm“®) means that no signifi¬ 
cant thermal or non-thermal X-ray emission is expected from it. By 
summing up (Pxray — Pradio)V from each region, where V is the 
volume of each region, we estimate that the total energy (Table 
required to heat all the entrained material in the lobe to 11 MeV is 
~ 9 X 10®’^ erg. 

Alternatively, pressure balance could be produced if the filling 
factor 0 < 1. This could occur in several ways, for example if we 
have overestimated the lobe volume or if significant variations in 
the density of relativistic particles and magnetic field exist within 
the lobes. 

The most likely location in which 0 may be < 1 is in regions 1 
and 5, which have values of (1 + A:)/0 much greater than those of 
the neighbouring regions. Since regions 1 and 5 of the lobes over¬ 
lap with the brighter X-ray emission north and south of the cool 
core, it is possible either that our assumption that the morphology 
of the lobe can be described as an annulus in this region is oversim- 
plistic, or that the lobe and IGM plasmas are mixed, with filaments 
or clouds of IGM gas penetrating the volume apparently occupied 
by the lobe. If we estimate the value of k based on the increase in k 
from regions 1-3 and 5-7, we find that we would expect k ~24 for 
region 1 and k ~23 for region 5. This suggests that 0 could be as 
low as ~0.59 for region 1 and ~0.48 for region 5. 

If thermal plasma is present in the lobes, we would expect it 
to cause Faraday rotation. We can therefore compare our expected 
density of entrained material with limits derived from the rotation 
measures for 3C 270 found by| Guidetti et al.| (|2011^ and |Laing & | 
|Bridle| ( [2014| l. While these measurements include the effects of ex¬ 
ternal Faraday rotation, they do provide an upper bound on the in¬ 
ternal rotation measure. Following |Feretti et al.| ( |199^ we calculate 
the expected rotation measure from: 
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Table 8 . Energy estimates for the particle content of the radio lobes 


Region 

Ee 

( 10®® erg) 

Rxray 

cm-2) 

Vtot 

( 10®^ cm®) 

(l + k)/4, 

Total energy 
( 10^® erg ) 

West Lobe 

1 

1.1 

30.2 

5.1 

42.5 

14.6 

2 

2.3 

11.8 

11.3 

17.4 

11.9 

3 

9.2 

9.1 

35.5 

10.5 

26.2 

4 

2.3 

4.8 

7.5 

A.l 

2.1 

East Lobe 

5 

0.6 

30.2 

3.1 

49.8 

9 

6 

1 

11.8 

4.8 

17 

5.1 

7 

7.1 

9.1 

28.7 

11.1 

21.4 

8 

5 

4.8 

13.7 

3.9 

3.3 


RM = 812 


/' 


neBz dz 


(4) 


where RM is the rotation measure in rad m“ , ne is the elec¬ 
tron number density in cm“^, Bz is the magnetic field along the 
line-of-sight in fj,G and L is the integration path (the depth of the 
lobe along the line-of-sight) in kpc. Based on the Guidetti et al. 
and Laing & Bridle results, an RM of order 10 rad m~^ seems 
representative for the system. We find that an electron gas den¬ 
sity of ~8.2xl0“® cm“^ would be required to produce RM=10 
rad m~^ purely from internal depolarization. The expected elec¬ 
tron density of thermal plasma entrained in the lobes (~1.4x 10“^ 
cm“®) is low compared to this limit, and would be expected to pro 


duce RM~0.02 rad m , well below the sensitivity of the Guidetti 


|et al.|p011^ maps. We therefore agree with the conclusion drawn 
by Guidetti et al., that the RM structures they observe are probably 
caused by foreground depolarization. 

We can also test whether regions 1 and 5, which we earlier 
suggested have filling factors less than unity, could contain a mix 
of thermal and relativistic plasmas, or whether it is more like that 
our adopted geometry leads us to overestimate their volume. Using 
eq uation and the IGM electron density derived from the X-rays 
in |0’ Sullivan et aLl ( |2011^ we find that we expect an internal RM 
of ~ 500 rad m“^ for both regions. While Guidetti et al. j201l) 
do find evidence of enhanced depolarization at the inner edges of 
the lobes, our predicted value is a factor ~20 higher than they ob¬ 
serve. This indicates that the relativistic and thermal plasmas along 
these lines of sight are probably not well mixed. This could suggest 
that the geometry we have assumed for the lobes close to the core 
is too simple, for example if filaments or sheets of thermal plasma 
penetrate the lobe while remaining magnetically separate. Alterna¬ 
tively, the internal properties of the lobe plasma may differ in these 
regions from those seen through the rest of the lobes. 

In |0’Sullivan et al.| ( |2011^ the total enthalpy, 4PV - the energy 
required to expand the lobe to its current size against the pressure of 
the IGM, plus the energy stored in the magnetic field and particles 
in the lobe - was estimated to be ~ 2.4 x 10®® erg. Our energy 
estimate of ~ 9 x 10®^ erg required to heat any entrained material 
is about 40% of the enthalpy estimated in |0’Sullivan et alT] ( |201 l| l. 
This is a factor of 1.8 greater than estimated by O’Sullivan et al., 
mainly owing to our more accurate treatment of lobe volume. The 
estimate changes little if the lower filling factors in regions 1 and 5 
are adopted. 

In summary, we have examined the constraints which can be 
placed on the particle populations in the radio lobes from the avail¬ 


able radio and X-ray measurements. We find that only a relatively 
small mass of entrained gas is needed to bring the lobes into pres¬ 
sure equilibrium with their surroundings, that this is probably con¬ 
sistent with entrainment of stellar material lost from stars within the 
jets, and that such material would have a low density and therefore 
would not have detectable effects on the polarization. 


7 CONCLUSIONS 

In this paper, we have presented a detailed spectral analysis of the 
FRI radio source 3C 270, emanating from the nearby group-central 
elliptical galaxy NGC 4261. We have analysed 240 MHz and 610 
MHz data from the GMRT, and data from the VLA at 1.55 GHz 
and 4.86 GHz. 

Combining the analysed GMRT and VLA data with data flux 
density measurements from the literature between 22 MHz and 
33 GHz, we used a Jaffe & Perola (IP) model to determine the 
asymptotic low-frequency spectral index of 3C 270 to be ainj = 
0-53^o.o 2- From the GMRT 240 MHz and VLA 1.55 GHz and 4.8 
GHz, we produced spectral index maps of the source which suggest 
a spectral index of a ~0.16 for the unresolved and probably self- 
absorbed core, a constant spectral index along the jets of a ~0.5 
and a gradual steepening from the tip of the jets through to the back 
of the lobes of a ~0.8—1. 

Assuming equipartition of energy between the magnetic field 
and the relativistic particles, and adopting a minimum Lorentz fac¬ 
tor of 7 min =100 for the radiating electrons, we estimate the radia¬ 
tive age of 3C 270. We modelled the observed 1.55-4.86 GHz spec¬ 
tral index trend in the lobes, fitting a Jaffe & Perola (IP) model. We 
find a radiative age of frad ~ 29 and 37 Myr for the west and east 
lobes respectively. This is about half the dynamical age of <75 
Myr found previously using X-ray data jo ’Sullivan et al.||20TT| ). 
We considered three possible scenarios for this discrepancy. The 
least likely is that the source underwent an early period of super¬ 
sonic expansion. While the shocks from such an expansion might 
now be outside the field of view of the available X-ray observa¬ 
tions, they would have very significantly heated the IGM, implying 
unphysically low pre-shock temperatures. 

A second alternative is that a backflow within the lobes may 
be transporting young, energetic electrons from the jet into regions 
of older plasma, leading us to underestimate the age. This can¬ 
not be ruled out, but previous high resolution observations of the 
base of the jets do at least argue that no backflow is observed on 
small scales. Finally, the third possibility, which we favour, is the 
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observed source is the product of multiple AGN outburst or large 
change in jet power, with the radiative age corresponding to the lat¬ 
est outburst or period of high jet power. [Worrall et ^ ( |2010[ l noted 
structures in the X-ray jet that could be indicate variations in power 
on 10"^ yr timescales, and larger scale variations in jet brightness are 
observed in the radio. A renewed outburst would not only inject a 
fresh population of young electrons into the lobes, but would likely 
drive shocks into the only population, causing reacceleration. The 
end result would be a source whose lobe expansion time might be 
similar to the X-ray age estimate while the radiative age represents 
the onset of more recent activity. 

Assuming that the pressures in the radio lobes equal that of 
the IGM we calculate that the ratio of energies of non-radiating to 
radiating particles, k, ranges from ~ 4 — 24 from the tip of the 
jets to the lobes close to the core (Table Fig.|10[>. This suggests 
that gas entrainment was more effective at earlier times, when the 
jets were shorter. Assuming, as a limiting case, that currently only 
material lost from stars in the jets is entrained, we estimate that the 
temperature of the entrained material would need to be 11 MeV if 
it is to bring the lobes into pressure balance with the IGM. This 
corresponds to a total energy of ~ 10®® erg required to heat the 
entrained material in the lobe, ~40% of the total enthalpy of the 
lobes estimated in |0’Sullivan et all]j2011| |. 

Based on our estimates of the density of the entrained mate¬ 
rial, we predict the likely internal rotation measure arising from this 
component, and compare it to RM measurements from the litera¬ 
ture. The density of ~ 1.4 x 10“^ cm~® indicates that the con¬ 
tribution of the entrained material to the Faraday depolarization is 
undetectable. In regions 1 and 5 where the thermal X-ray emitting 
plasma co-exists with the relativistic plasma, we would expect a 
RM of ~ 500 rad m~^ for both regions if mixing had taken place. 
This is ruled out by the observations, implying that little mixing 
between the thermal and relativistic plasmas has occurred. 
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